Matter-wave bistability in coupled atom-molecule quantum gases 
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We study the matter-wave bistability in coupled atom-molecule quantum gases, in which het- 
eronuclear molecules are created via an interspecies Feshbach resonance involving either two-species 
Bose or two-species Fermi atoms at zero temperature. We show that the resonant two-channel Bose 
model is equivalent to the nondegenerate parametric down-conversion in quantum optics, while the 
corresponding Fermi model can be mapped to a quantum optics model that describes a single-mode 
laser field interacting with an ensemble of inhomogeneously broadened two-level atoms. Using these 
analogy and the fact that both models are subject to the Kerr nonlinearity due to the two-body 
s-wave collisions, we show that under proper conditions, the population in the molecular state in 
both models can be made to change with the Feshbach detuning in a bistable fashion. 
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I. INTRODUCTION 

We study the matter-wave bistability in coupled 
atom-molecule quantum gases, in which heteronuclear 
molecules are created via an interspecies Feshbach res- 
onance involving either two-species Bose or two-species 
Fermi atoms at zero temperature. We show that the res- 
onant two-channel Bose model is equivalent to the nonde- 
generate parametric down-conversion in quantum optics, 
while the corresponding Fermi model can be mapped to a 
quantum optics model that describes a single-mode laser 
field interacting with an ensemble of inhomogeneously 
broadened two-level atoms. Using these analogies and 
the fact that both models are subject to the Kerr nonlin- 
earity due to the two-body s-wave collisions, we show that 
under proper conditions, the population in the molecu- 
lar state in both models can be made to change with the 
Feshbach detuning in a bistable fashion. 

The ability to cool and trap neutral atoms down to 
quantum degenerate regime has created a host of new and 
exciting problems that are increasingly interdisciplinary, 
bridging in particular the atomic, molecular, and opti- 
cal physics and the condensed matter physics. The rich 
knowledge and experience accumulated over the past sev- 
eral decades in these fields have dramatically accelerated 
the progress of ultracold atomic physics. An example 
that serves to illustrate how the interdisciplinary fields 
learn and benefit from each other is the phcnomonon of 
atomic pairing where a bosoinc molecule is coupled to 
two bosonic or fermionic constituent atoms via Feshbach 
resonance or photoassociation. So far this is the only vi- 
able approach to create ultracold molecules. It is also 
an ideal test ground for studying coupled atom-molecule 
condensates and the BCS-BEC crossover [2j. The lat- 
ter is thought to be underlying the mechanism of high 
temperature superconductors and extensively studied in 
the realm of condensed matter physics. In addition, the 
coupled atom-molecule systems have deep quantum opti- 



cal analogies [1, [B[ : bosoinc molecules coupled to bosonic 
atoms (which we will refer to as the bosonic model in this 
paper) is the matter-wave analog of parametric coupling 
of photons which has important applications in generat- 
ing nonclassical light fields and, more recently, in quan- 
tum information science; while the system of bosonic 
molecules coupled to fermionic atoms (which we will refer 
to as the fermionic model) can be mapped to the Dicke 
model where a light field interacts with an ensemble of 
two-level atoms, a model having fundamental importance 
in the field of quantum optics. 

In this work, we will further explore these quantum 
optical analogies of the atom-molecule system and fo- 
cus on the important effects of binary collisional interac- 
tions between atoms which are largely ignored in previous 
studies 0, [j| . We show that the atom-atom interaction 
introduces extra nonlinear terms which, under certain 
conditions, give rise to matter- wave bistability in both 
bosonic and fermionic models. Hence, we may establish 
the connection between the coupled atom-molecule quan- 
tum gases and the nonlinear bistable systems Q that 
have been extensively studied in the 80 's in the context 
of nonlinear optics, due both to its fundamental inter- 
est, and to its many practical applications in fast optical 
switches, optical memory, laser pulse shaping, etc. 



II. BOSONIC MODEL 

In what we call the bosonic model, a molecule as- 
sociated with annihilation operator a m is coupled to 
two non-identical atoms labeled as | f) and | J.) with 
corresponding annihilation operators cif and a^, respec- 
tively. Here we consider two types of atoms in order 
to make direct comparisons with the fermionic model 
to be treated in the next section, for which only unlike 
fermionic atoms can pair with each other and form a 
bosonic molecule. Futhermore, in this work we only con- 
sider zero-temperature homogeneous case so that all the 



2 



bosons are condensed into zero center-of-mass momen- 
tum states. 

The second quantized Hamiltonian reads 

H = 5 a) m a m + g (a] n a ] a i + h.c.) + x«aJ S i 6 .A > C 1 ) 

where the detuning 8 represents the energy difference 
between the molecular and atomic levels which can be 
tuned by external field, g is the atom-molecule coupling 
strength and Xij = Xji is the s-wave collisional strength 
between modes i and j. This system has been studied in 
Ref. 0]. For completeness and better comparison with 
the fermionic model, we briefly state some of the main 
results relevant to the focus of this work — matter-wave 
bistability — and direct readers to Ref. Q for more de- 
tails. 

For our purpose, we take the standard mean-field ap- 
proximation and replace operators dj with c-numbers 



"] — V Nj e lVj . The mean- field Hamiltonian takes the 
form: 

H = 2k{y 2 -y) + 2vy+{\-2y)^hj cos^, (2) 

where 

y = 0.5 [1 - (JVf + iVj /N] = N m /N , cp = ^+ Vl -<p m , 

are a pair of conjugate variables, representing the molecu- 
lar population and phase mismatch, respectively. Other 
quantities are defined as 



G = gV2N , 

A = N (XTT + XU + Xrnm + 2\ U 

v = [6 + XTT + XU + (N -1) Xm: 



2XmT - 2 Xml)/G, 

- AT XmT - Nxmi\ /G , 



with TV = + N± + 2N m a constant of motion represent- 
ing the total number of atoms, and we have assumed that 
the number of atoms in states | f) and | J.) are equal, i.e., 
iVf = N±. In addition, we will focus on the stationary 
states with <p = 7r which has lower energies than the ones 
with (p = 0. 



A. Quantum Optical Analogy 

It is quite clear from the form of the second-quantized 
Hamiltonian in Eq. ([T]) that without the collisional terms 
our model will reduce to the trilinear Hamiltonian de- 
scribing the nondegenerate parametric down-conversion 
in quantum optics [H, Q. In this analogy, the molecu- 
lar mode plays the role of the pump photon, where the 
two atomic modes are the signal and idler photons, re- 
spectively. The collisional terms would correspond to the 
Kerr-type cubic nonlinearity which will be present in the 
optical system if the light fields propagate in some non- 
linear medium [10(. 



B. Bistability 

In the absence of the collisions or Kerr nonlinearity 
(i.e., A = 0), the system does not exhibit bistability. 
This can be seen by studying the properties of the mean- 
field Hamiltonian H in Eq. which can be simplified 
as (taking ip = n) 

H = 2vy-(l-2y)^Zy~, (3) 

The stationary state correspond to the solution of 

dH 



dy 



= 2^ + 3^-1/^=0. 



(4) 



For a given detuning is, the stationary state is unique: 
0.5, u<-l 



U-v + ^TZf, V >-1 



(•5) 
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For A 7^ 0, using Eq. ([2|), the stationary condition is given 
by 



dH 
dy 



2i/ + 3^/21/ -l/V2j/ = 0, 



where we have defined 



v' = v + k(2y- 1) 



(6) 



(7) 



Note that Eqs. (|3|) and ^ have the same form. In other 
words,we can express the effect of collisions as a nonlinear 
phase shift for molecules that modifies the detuning v. 
Consequently, the stationary solution for A ^ should 
have the same form as in Eq. ([5]) but with v replaced by 
v' , which makes yo an implicit function of the detuning 
v. To find the explicit dependence of yo on v, we can 
use the graphic method as illustrated in Fig. [1] For the 
example given, we obtain three stationary states. Further 
analysis shows that the middle solution is dynamically 
unstable and the other two are stable solutions 7] . Such 
a behavior is typical in bistable systems Q . 

The graphics of Fig. [1] also shows that, in order to have 
multiple stationary solutions, the slope of the straight 
line (given by 1/2A) must be negative and cannot be too 
steep. More specifically, the slope of the straight line has 
to be larger than the slope of the curve at v — —1, and 
this leads to the condition 

A<-1, (8) 

in order for the system to exhibit bistability. 



III. FERMIONIC MODEL 

In the fermionic model, we denote as the anni- 
hilation operator for an atom with spin a (=| , J,), mo- 
mentum fik, and energy = h 2 k 2 /(2m), and as before 
denote the annihilation operator for a molecule in 
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FIG. 1: For given A and u, the thick solid line represents 
ya(v') and the thin dashed straight line represents Eq. (0. 
The intersects are the stationary solutions. Here we take 
1/2A = -0.1 and v = 0.4A. 

state | to) with zero momentum, the second quantized 
Hamiltonian reads: 

H = Y1 e fc»L«k CT + U «kT a -k+qi«-k'+qiak'T 

k.tr k,k',q 

+ va\ n a m + —j= Y (aLa-kiflkT + h.c.) , (9) 

where V is the quantization volume. Hamiltonian §§§ has 
the form of the two-channel model of BCS-BEC crossover 
where only the condensed molecule part is considered 
fill ] . Following the Hartree-Fock-Bogoliubov mean- field 
approach by dividing the two-body collision into a 
part related to the BCS gap potential A — Up, where 

P = X! («-kJ.akT) /V, 

k 

and a part related to the Hartree potential 

V h = Uj2(aL^*)/( 2 V), (10) 

k<r 

where we again assume equal population in | f) and | J.) 

atomic states, i.e., (fljL-OkT) = (a^aki), we may express 
the Hamiltonian as 

k,cr 

+ [( Up + 9a m /^/v) a k T«-ki + h - c ■ ( n ) 

k 

Defining N = 2a\ n a m + ^ k a a^a^ as the number op- 
erator which is a constant of motion, we may rewrite the 
term proportional to Vh in (fTTjl as 

= V h N - (Un-2U{tfb)/vYtfb, (12) 



where n — (N) jV is the constant total atom number 
density. In our derivation, Vh arises from the two-body 
atom-atom collision. In general, additional terms repre- 
senting atom-molecule and molecule-molecule collisions 
are also present. These additional terms will modify the 
coefficient U in the definition of Vh [Eq. flit)))], which is the 
counterpart of A in the bosonic model, but the general 
form of Eq. (|12p will remain valid. In the following, we 
will refer to this term as the collisional term. Through 
Eq. (fT2|) . we have expressed the effect of the two-body 
collisions as a nonlinear energy shift of the molecules 
(along with a constant energy bias VhN), in complete 
analogy with the bosonic model. We remark that in the 
usual one-channel model of the mean-field BCS theory 
valid when the molecular population is negligible, the 
collisional term just represents an unimportant constant 
energy shift. 

As usual, dko- (t) and a m (t) obey the Heisenberg equa- 
tions of motion based on Hamiltonian (|11| . By replacing 
Bose operator a m with the related c-number c — (b) / \fv 
and Fermi operators Ska- (t) with the familiar Uk (t) and 
Vk (t) parameters through the Bogoliubov transformation 

Okt = u* k dkt + Vk a [ _ kl and aL ki = -v* k a k | + u k aL ki , 
where dko- are the Fermi quasiparticle operators, we ar- 
rive at the following set of mean-field equations of motion 

ihc = v e c + gp, (13a) 
iMk = -tkUk + AgWfe , (13b) 
ihvk = A e iik + tkVk , (13c) 

where p = u%v k /V, A e = gc + Up , and 

v e = v- Un + 2U\c\ 2 , (14) 

is the effective detuning which contains a Kerr nonlinear 
term 2U\c\ 2 whose origin can be traced to the two-body 
collisional shift. This set of equations describes the dy- 
namics at zero temperature where the state of the system 
can be described as the quasiparticle vacuum. 

A. Quantum Optical Analog 

In several previous studies where the collisional term 
is neglected, it has been pointed out that the fermionic 
model can be mapped to the Dicke model in quantum 
optics [1, [13] as schematically shown in Fig. [5] (see below 
for details). In fact, this model was recently shown to 
display collective dynamics similar to photon echo and 
soliton-like oscillations in transient collective coherent 
optics [H[ . Such a connection can be traced to the work 
of Anderson's spin analogy [l4| for the BCS problem. 

To show what is the quantum optical analogy of the 
collisional term, let us rewrite Eqs. (|13)) in a form more 
familiar in cavity optics. To this end, we first introduce 
a set of new variables 

Pk = 2u* k v k , D k = \u k \ 2 - \v k \ 2 , £ L = 2iA e , 
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FIG. 2: (Color online) Mapping of the two-channel resonant 
Fermi superfluid model to the Dicke model. The bosonic 
molecules and the fermionic atoms in the former are mapped 
to the cavity laser field and an ensemble of two-level atoms in 
the latter, respectively. See text for details. 

and recast Eqs. (|13b| . (|13c|) into 

hP k = -i2e k P k - £ L D k , (15a) 
hD k = (£* L P k +£ L P*)/2. (15b) 

Interpreting P k and D k as the microscopic polarization 
and population inversion, respectively, Eqs. (JTSJ) then be- 
come the optical Bloch equation that describes the in- 
teraction between a local electromagnetic field £l and 
a fictitious two-level atom, characterized with a transi- 
tion energy 2e k [16]. This analogy is consistent with the 
fact that there exists a one-to-one mapping between pairs 
of fermion operators and Pauli matrices when the BCS 
pairing mechanism is taken into account [14| . 

In this optical analogy, the local electric field £l = 
£+£i contains two contributions because of A e = gc+Up. 
The first of these (£ — i2gc) is equivalent to an av- 
erage macroscopic field, whose dynamics is described 
by Eq. (|13a|) . which can now be interpreted as the 
Maxwell's equation for the cavity field £ with cavity 
detuning v ei driven by a macroscopic polarization den- 
sity p = J2k-Pk/(2V) of an inhomogeneously-broadened 
medium [see Fig. [SJ. The second part £i = i2Up may be 
regarded as the internal field at the atom due to the col- 
lective dipole polarization of the nearby two-level atoms 
in the ensemble. As such, £l = £ + £i here bears a 
direct analogy to the Lorentz-Lorenz relation in optics 
[l7l ]. Note that had the collisional term been neglected 
(i.e., U = 0), there would have been no internal field 
contribution, nor would there have been the Kerr nonlin- 
earity in the equation for the bosonic mode. For U 7^ 0, 
both of these terms will be present. Under such a cir- 
cumstance, Eqs. (|13al) and (fTS")) represent the generalized 



optical-Bloch equations in which the Lorentz-Lorenz re- 
lation is explicitly incorporated [18j, and hence can lead 
to interesting nonlinear phenomena just as they do in 
optical systems. 



B. Bistability 

Having established this analogy, we now look for the 
steady state solution from Eqs. (|13ap and (fT5|) . As is 
well-known, the operation frequency of a laser field is 
not known a priori; but is established through the so- 
called mode pulling — the dynamical competition be- 
tween atomic and cavity resonances. A similar argument 
holds for the molecular field c. For this reason, we adopt 
the following steady-state ansatz 

c _ ce -«/**/A Pk Pk e -*i*/* t Dk - D k 

where the same symbols are used for both dynamical 
and steady-state variables for notational simplicity. The 
molecular chemical potential, 2[i, is just the correspond- 
ing lasing frequency in the cavity optics model. From the 
steady state equations obtained by inserting this station- 
ary ansatz into Eqs. ()13a[) and (fT5"]) . we can easily find 
that (a) there always exists a trivial solution or a "non- 
lasing" state with A e = or equivalently c = 0, which 
corresponds to the non-superfluid normal Fermi sea; and 
(b) a non-trivial solution with its /x, A e and c determined 
sclf-consistently from the gap equation 

1 1^1 

U - g*/(ve - 2/i) ~~W^E~ k ' {b) 

with E k = \J (e k — /i) 2 + A 2 , the number equation 

and an auxiliary relation 

| 5 A e | = \c{u e - 2y)[U - g 2 l{v e - 2 M )]| . (18) 

The integral in the gap equation (fT6|) under the assump- 
tion of contact interaction is known to be ultraviolet di- 
vergent. To eliminate this problem, we renormalize the 
interaction strength U and g, as well as the detuning v 
in (|T6|) . while U in the collisional term is rep laced by the 
background interaction strength Uq [HI, \2(\ . 

Note that there exists, in the single-mode inhomo- 
geneously broadened laser theory [21(, a similar set of 
steady-state integral equations, which, due to lasers be- 
ing open systems, are obtained under different consider- 
ations. For example, the requirement that the cavity loss 
balance the saturated gain leads to the "gap" equation, 
whose primary role is to limit the laser intensity; while 
the phase matching condition translates into the "num- 
ber" equation, whose main responsibility is to assign the 
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relevant equations of motion maintain the same forms if 
we substitute v by 



FIG. 3: (Color online) Free energy density / as a function of 
A e and |c| 2 at v = 0.02 (a) and v — 0.2 (b). Extremum points 
are indicated by 'x' (minimum) and '+' (saddle point). /, 
A e , and v are all in units of Ef = (3n 2 n) 2 ^ 3 / (2m), the Fermi 
energy of the non-interaction system. In all the examples 
shown in this paper, the physical parameters corresponding 
to go and Uo are 1.2 Ef [k 3 p and —60.7 Ep/k F , respectively. 
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FIG. 4: Molecular population |c| 2 as a function of detuning. 
Vertical line in (a) indicate the critical point of a first-order 
phase transition. In (a) the collisional term is included while 
it is neglected in (b). 



amount of mode pulling of the laser field relative to the 
cavity resonance. 

An alternative way to derive Eqs. (fT!)|) - (fT5|) is from 
the energy density. The zero-temperature energy density 
/(A e ,c, n) = (H)/V can be calculated using Hamilto- 
nian (fTTj) and the Bogoliubov transformation as [2(| 



e k - a - E k (A e - gc) 2 



V 



U 



+ (y e - 2/i)|c| 2 + fm. 



(19) 

The extremum conditions df/dA e = df/dc = 0, lead 
to Eqs. (flB]) and (fT5|) , respectively, while the condition 
df /da = results in the number equation (jTTJ). 

Figure [3] illustrates the energy density in the |c| 2 -A e 
plane for different detuning v. For any given pair of (c, 
A e ), fi is calculated self-consistently using the number 
equation (|17p . Typically, / has only one extremum which 
is a minimum point as shown in Fig. [3£a). However, 
in the regime v E (—0.08, 0.13) Ep, / possesses three 
extrema: two of them are local minima and the third 
a saddle point. An example with v = 0.02 is shown in 
Fig. [Ib). 

To gain more insights into the bistable behavior, we 
may carry an analogous analysis as in Sec. Ill Bl In the ab- 
sence of the collisional term, steady-state molecular pop- 
ulation |c| 2 is a smooth monotonically decreasing func- 
tion of v and the system does not exhibit bistability: 
As v increases, molecules decompose into atoms. This is 
shown in Fig.2£a). When collisional term is included, the 



2U \c[ 



(20) 



Hence the solution |c| 2 as a function of v' is represented 
by the same curve as in Fig. 0Ja). To find |c| 2 as a func- 
tion of v, we need to find the intersects between this curve 
and the straight line representing Eq. (j2"U)) . In direct anal- 
ogy to the graphic method in Fig. [TJ for Uo sufficiently 
large and negative, these two curves have three inter- 
sects and the system exhibits bistability. One example is 
shown in Fig.fJJb). The vertical line in Fig.QJb) indicate 
the critical point of a first-order phase transition: across 
this line, the ground state jumps from the upper branch 
to the lower one. For the parameters used, this occurs at 
v c = -0.01 E F . 

To check the stability of these steady states, we have 
solved the dynamical equations (fll5|) using the slightly 
perturbed steady state solution as the initial condition. 
From the dynamical evolution of the system one can 
see that, just like in the bosonic model, the states in 
the upper and lower branches are dynamically stable: 
when slightly perturbed, they exhibit damped oscilla- 
tions around their equilibrium values. These oscillations 
can be further understood from the excitation spectrum 
of the corresponding steady state. This can be done us- 
ing a linear stability analysis, which is also the standard 
tool for studying laser instabilities [2l],|2^|. The spectrum 
is found to contain a discrete part which determines the 
oscillation frequencies, and a continuous part which con- 
tributes to the damping of these oscillations at a much 
longer time scale [2J|. By contrast, the states in the mid- 
dle branch are unstable as small perturbations will lead 
to large departures. 



C. Dynamics 

The bistability has important ramifications in atom- 
molecule conversion dynamics. When the collisional 
term is unimportant and negligible, one can easily create 
bosonic molecules from fermionic atoms by adiabatically 
sweeping the Feshbach detuning across the resonance. 
As long as the sweeping speed is sufficiently slow, the 
molecular population will follow the steady-state curve as 
shown in Fig.|4ja). By contrast, when bistability induced 
by the collisional term occurs, the adiabaticity condition 
will necessarily break down. Fig. O displays the dynami- 
cal evolution of the bosonic population when the detun- 
ing is swept starting either from a large positive or a large 
negative value. We can see that the steady-state curve 
can be followed up to the point where the stable states of 
the upper and lower branches and the unstable states of 
the middle branch join each other (indicated by v\ and v<i 
in Fig. [5]), where the population suddenly jumps between 
the two stable branches. Note that the critical detuning 
v c for the first-order phase transition as indicated by the 
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FIG. 5: (Color online) Dynamics of atom-molecule conversion 
as illustrated by the molecular population when the detuning 
v is slowly swept. Curve (a) is obtained by sweeping v from 
positive to negative values, while curve (b) is obtained by 
sweeping v in the opposite direction. The dotted line is the 
steady-state molecular population, the same as in Fig. |4jb). 

vertical line in Fig.[5Jb) lies between v\ and V2- The dy- 
namical population curve thus exhibits hysteresis in the 
vicinity of the first-order phase transition. In this way, 
by tuning the detuning in the vicinity of V\ or 2/2 > a- n 
atom-molecule switch can be realized. Similar behavior 
is also found in the bosonic model. 



IV. CONCLUSION 

In conclusion, we have studied the matter-wave bista- 
bility in coupled atom-molecule quantum gases in both 
the bosonic and the fermionic models. These two cases 
can be mapped to two very different quantum optical 
models: parametric downconversion in the former and 
generalized Dicke model in the latter. Nevertheless, one 
important common feature for both cases is that bista- 
bility can be induced by collisional interactions which 
give rise to Kerr nonlinearity. We hope that our work 
will motivate experimental efforts in demonstrating the 
matter-wave bistability we predicted here. 
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